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ABSTRACT 


The dispersion curves have been computed for a wide 
range of the wave length from the electromagnetic waves 
to the electrostatic waves in the case of the magnetoactive 
warm plasma with a Maxwellian velocity distribution func- 
tion. The computation was carried out mainly for the 
perpendicular propagation mode. The upper hybrid resonance 
is the connection point of the electrostatic waves and the 
electromagnetic waves. The electrostatic waves not associa- 
ted with the upper hybrid resonance are subjected to elec- 
tron cyclotron damping when the wave length becomes long. 
Oblique propagation is allowed for the electrostatic waves 
in a frequency range from the plasma frequency to the upper 
hybrid resonance frequency in the long wave length region 
where Landau damping can be neglected and where the electro- 
static mode smoothly connects to the electromagnetic X-mode . 
In a slightly inhomogeneous plasma the Bernstein mode elec- 
trostatic wave can escape by being converted into the O-mode 
electromagnetic wave; two reflections take place during this 
escape process. The frequency range of the escape coincides 
with the Cerenkov radiation of the electromagnetic waves. 


but these two mechanisms indicate clear contrast in its origin. 
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but these two mechanisms indicate clear contrast in its origin. 



The electromagnetic Cerenkov radiation is produced interact- 
ing with the high energetic plasma beam with velocity range 

100 v . , < v < 1000 v (for the nonrelativisitic plasma), 
th tjti 

where v is the electron thermal velocity, while the origin 

vll 

of the electrostatic waves is in the interaction with the 


low energetic plasma beam. 



1 . INTRODUCTION 


The Bernstein mode wave (Bernstein, 1958; Crawford, 

1965) is the electrostatic approximation of the dispersion 
equation of plasma waves in a magnetoactive warm plasma 
for perpendicular propagation. In the extreme case of 
long wave lengths, the plasma waves exhibit the well known 
electromagnetic characteristics which were originally in- 
vestigated by Appleton and Hartree (for a review, see 
Budden, 1961) for the case of a cold plasma. Shkarofsky 
(1968) and Dougherty (1969) have discussed the coupling 
region between the Bernstein waves and the electromagnetic 
waves in frequency ranges near the electron cyclotron 
harmonic frequencies. 

In the present paper, the dispersion equation is investi- 
gated numerically to clarify the full nature of the plasma 
waves from the electrostatic to the electromagnetic modes 
over a wide frequency range. An escape mechanism of the 
electrostatic wave as a form of the electromagnetic wave is 
proposed relating to the upper hybrid frequency. The detailed 
discussion is carried out relating to the escape mechanism 
when the wave packet propagates through the slightly in- 
homogeneous plasma. The connection of the dispersion relation 
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for the oblique propagation in a wave length which is long 
enough to be free from the Landau damping effect has also 
been investigated relating to the transfer of the wave 
energy from the X-mode to the O-mode wave. 

The extraordinary mode electromagnetic wave at the 
upper hybrid frequency has been observed in the laboratory 
plasma using a positive column of a weakly ionized hot- 
cathode arc discharge in helium with the magnetic field 
parallel to the column (Bekefi et al . , 1962). The reception 
of the electromagnetic wave which propagates perpendicular 
to the magnetic field strongly suggests that the electrostatic 
oscillation can escape at the upper hybrid frequency. The 
present work provides a theoretical base for interpreting this 
evidence in terms of the escape mechanisms of the electro- 
static wave as a form of the electromagnetic wave. 

In the last part of this paper a brief discussion con- 
cerning the discrimination between the electromagnetic 
Cerenkov radiation and the electromagnetic wave which is pro- 
duced as a reformation of the electrostatic wave propagating 
through an inhomogeneous medium is discussed because these 
two kinds of electromagnetic waves occur in the same 
frequency range between the plasma frequency and the upper 
hybrid frequency. 
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2. DISPERSION RELATIONS 

The dispersion relation which has been summarized by 
Stix (1962) for the magnetoactive warm plasma can be derived 
from the wave equation 


n x n x E - (K) E ~ 0 (1) 

where n = (c/iu) k, k is the wave vector, c is the velocity 
of light, and uj is the angular frequency. The matrix (K) is 
the equivalent dielectric tensor for the magnetoactive warm 
plasma with a Maxwellian velocity distribution function. 

The elements of the tensor (K) have the form 
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and j indicates the different species of particle in the 

plasma. In the above expression e . = -1 for electrons and 

J 

+1 for ions; also rr . , 0., T , m., k , k , and v are the plasma 

frequency, the cyclotron frequency, the temperature corresponding 

to the velocity component perpendicular to the magnetic field, 

particle mass, the z and x components of the wave vector k, 

and the z component of the particle velocity v, respectively. 

The function X in eq. (2) is the modified Bessel function 

n 

with argument (=k «T/m . ) and I’ represents its first 

j x j n 

derivative. In the case of T = T = T and zero drift velocity 

i. a 

the function in eq. (2) , for a sufficiently high frequency 
range where only the electron motion is important, is re- 
written by 
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where m, and ? are the electron mass and a real variable, 


respectively. In the present case of nearly perpendicular 
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propagation « n » 1; hence we can use the following approximate 


expression for F q : 
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with an error less than 0.2% for > 3 . 
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The wave equation in the source free region, i.e. eq. (1), 
can be rewritten as 
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For the case of the nearly perpendicular propagation, 

n - 0 , the wave can only exist for the condition, 
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This relation Is rewritten in a more simple form as 
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2 

where A ~ (c m/WT) and q = a )/Q. The values of A in the ionosphere 
(T 2000°K) , in the magnetosphere (T = 10^ °K) , and in the solar 
corona (T = 10® °K) are 2.967 X 10® , 5.94 x 10 4 and 5.94 x 10 2 3 , 
respectively. Thus, the only solution of eq. (7) in the 
short wave length domain, i.e., large x, is 

K - 0 (k) 

xx 


6 


This relation is the dispersion equation of the well known 
Bernstein mode electrostatic waves. 

In the extreme case of the long wave length, the elements 
of the dielectric tensor become 



- 1 ), 


where 



(rr/O) 2 


(9) 


Rewriting k c/uj by n, eq- (7) reduces to the 

X 

Appleton-Hartree relation for 0 = 90° 

2 2 , , 2 
n = 1 - q N /(q -1) 

and (10) 

2 2,2 
» = 1 - q N /a . 

In the case of intermediate wavelengths all of the terms in 

eq. (7) must be retained. The next section presents numerical 

solutions of eq. (7) corresponding to this case. 

3 • N UMERICAL SOLUTIONS OF THE FULL WAVE DISPERSION RELATION 

— _ 6 

Computations have been carried out for A - 2.97 x 10 

(T ~ 2000°K) and several values of the plasma parameter 

q (= tt/H) ranging from 0.2 to 4.5, as indicated in Figs. 1 
N 

a»3 

to 4. The curves with k R < 5 x 10 are identical with the 

x 
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dispersion curves of the electromagnetic waves except for 

the frequency range close to the cyclotron harmonics. The 

-1 

curves with k R > 10 are identical with the dispersion 
x 

curves of the Bernstein mode electrostatic waves where 



In the intermediate wavelength region, i.e., between 
the extremes of pure electrostatic waves and pure electro- 
magnetic waves, only waves described by the dispersion curves 
near the upper-hybrid frequency can persists; the waves 
corresponding to the other dispersion curves, which merge 
into the harmonics of the electron cyclotron frequency, are 

absorbed due to strong cyclotron damping. Although k = 0 

z 

is assumed in the present investigation, the actual propaga- 
tion modes will include infinitely small k^ values rather 
than zero. This means that the wave is subjected to strong 
cyclotron damping effects when < 3 (see F q in eq. (3)) 

which corresponds to the frequency range. 

nfi - 3k < to < nn + 3k J^*2. . 

z i m z i m 

Near the upper hybrid resonance, however, the wave is undamped 
and it is indicated (see Figs. 1 to 4) that the hybrid frequency 

is neither the resonance (k -* *», i.e., the wave length tends 

to zero) nor the cut off (k 0, i.e., the wave length tends 

to infinity), but provides a wave which creeps with relatively 


( 11 ) 
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low group velocity (dw/dk ), When q < 0.8, the minimum 

X JN 

value for the group velocity at the hybrid resonance frequency 

is less than 0.04 v iL where v.. is the electron thermal velocity 

th th 

which is defined by RQ. When q„ > 1.6, the minimum value for 

N 

the group velocity is between 0.08 and O.lv^. The magnetic 
field component gradually increases for wavelengths greater 
than about 200 nR which can be considered as the transition 
region between the electrostatic wave and the electromagnetic 
wave . 

4. ESCAPE MECHANISMS OF THE ELECTROMAGNETIC WAVES FROM THE 
THERMAL PLASMA 

a) Starting point 

The above evidence that the electrostatic wave near 
the upper hybrid resonance continuously changes into an elec- 
tromagnetic wave as the wave length increases, provides a 
theoretical foundation for the escape mechanism of Bernstein 
mode electrostatic wave energy as electromagnetic wave energy . 
These electrostatic waves may originate from instabilities in 
the plasma. In order to obtain the variation of the wave 
length, it is necessary to have an inhomogeneous plasma. To 
illustrate the escape mechanism we will consider a condition 
in which the direction of the electron density gradient makes 
a finite angle with respect to the superimposed magnetic field, 
and the wave starts propagating in the direction perpendicular 
to B corresponding to decreasing q N (see Fig. 5). Let ns 
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consider a test wave packet that is constructed by a 
narrow band of sinsoidal waves with a center frequency uu^. 

These waves are indicated on the dispersion curve in Fig. 

6 by A^. Assume that this wave packet starts from the 
point A in the inhomogeneous plasma depicted in Fig. 5. 

When the wave is electrostatic, the propagation direction 
is strongly restricted by the direction of the constant 
magnetic field; since the component parallel to the magnetic 
field is subjected to thermal plasma damping the propagation 
direction is strongly localized in the direction perpendicular 
to the magnetic field. 

b) The turn of electrostatic wave 

Since we are assuming that the wave initially propagates 
in the direction of decreasing q._, the dispersion curves will 
flatten out (see Figs 1-4) and the point on the curve corres- 
ponding to the center frequency of the wave packet will 
eventually merge with the point in Fig. 6. 

At this point the group velocity of the wave change 

from Bu)/5k < 0 to duj/dk > 0 and the test wave packet is re- 
x x 

fleeted back towards the high q N region. This turning point 

is shown by the point B in Fig. 5. 

m 

c) Passing through the upper hybrid resonance 

After turning the wave packet propagates toward the 
region of the higher q N where the dispersion curves build up 
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(see Figs 1-4) and the point on the curve correspond to 
the center frequency of the wave packet moves toward in 
Fig. 6, where marks the level of the upper hybrid resonance 
frequency. The group velocity of the wave around the upper 
hybrid resonance is not zero and the wave packet keeps pro- 
gressing toward the high q^ region. When the wave length 
becomes larger than 200 ttR (for example), the organized motion 
of the plasma generates a current which produces a significant 
magnetic field. Hence, the wave continuously changes from 
electrostatic to electromagnetic as the wave length increases 
while the packet propagates toward the high q^ region. For 
this case, the corresponding points are indicated as and 

C in the model plasma in Fig. 5, and on the dispersion curves 
d 

in Fig. 6, respectively. 

d) Angle change and mode conversion 

After passing through the region of the upper hybrid 
resonance frequency the ray direction of the test wave packet 
becomes more or less density dependent. This is partly because 
the wave propagation in an oblique direction with respect to 
the magnetic field is not severely damped in the long wave 
length region. In Fig. 7, the dispersion curves are plotted 
for oblique propagation with angles 0 = 0°, 15°, 30°, 45°, 60°, 
75° and 90° for q = 1.6. For oblique propagation, plasma 
thermal damping is very severe in the short wave length (large 
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k value) region. Thus the curves are plotted only for wave 

X 

lengths longer than 26.7 {R cos 0/(q + n) } . This criterion 

of the wave length limit is chosen to make a less than 3, 

n 

in eq. (3), which makes the imaginary part of the function 

-4 

i of F less than 5 x 10 . The value k for a given constant 

no x 

frequency becomes smaller and the group velocity becomes 

larger, with increasing q , until the center frequency of the 

wave packet coincides with (in Fig. 6). Then, the path 

of the wave packet propagation starts bending as indicated 

schematically in Fig. 5, after passing through the point C . 

When the wave packet arrives at the point D where a> = n, 

mo 

all the dispersion curves for the different angles merge 

together as indicated by the point in Fig. 6. At this 

point the part of the wave energy that can be transformed into 

the O-mode wave escapes from the plasma region from the point 

D as indicated in Fig. 5; the other part of the wave energy 
m 

propagates in the plasma until it is reflected from the cut- 
off frequency u> for the Z wave. After reflection the above 

Ju 

mentioned mode conversion again takes place when the point 
D d of Fig. 6 is encountered (E^ of Fig. 5). Part of this 
wave energy escapes as an O-mode wave and part of it experiences 
the reciprocal process of mode conversion, i.e., electro- 
magnetic to electrostatic. This reciprocal process corresponds 
to the absorption of electromagnetic waves by a plasma (Budden 
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(1961), Denisov (1958) and Stix (I960)). 

RELATION WITH CERENKOV RADIATION 

As discussed in the previous section, electrostatic 
wave-energy can escape from an inhomogeneous plasma in the 
form of electromagnetic energy in the region where the wave 
frequency is between the local plasma frequency and the 
local upper hybrid frequency. This frequency range coincides 
with the frequency range of observed electromagnetic wave 
Cerenkov radiation which has been investigated by several 
workers (Jelly (1958), Kolomenskii (1956), Eidman (1962), 
Mckenzie (1964) , Liemohn (1965) , Seshadri and Tuan (1965) , 
and Bauer and Stone (1968)). 

The condition for Cerenkov radiation is 

v cos 0 = u>/k 

where v, and 0 are the scaler velocity of the electron beam 

and the angle between the propagation vector k and the elec- 

tron beam velocity v, respectively. The above mentioned 
— > 

authors let k -* 00 for the non relativistic beam where 
v /c o. But as has been confirmed in the present paper, 
the electromagnetic approximation to the dispersion equation 
breaks down when k -* 00 . A proper solution requires the 
consideration of the full dispersion equation together with 
eq. (12) . This equation is rewritten as 


( 12 ) 
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k R 
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q 

(v/v., ) cos 9 
th 

The solution of eq. (13) can be obtained graphically as 

indicated in Fig. 8; the diagram indicates that Cerenkov 

radiation occurs for different wave lengths corresponding 

to different cross points of a given (v/v, . ) cos 9 curve 

th 

and the dispersion curves. 

From this process we can see that the Cerenkov radia- 
tion associated with a high speed beam, (v/v ) cos 0 > 100, 

L XI 

is electromagnetic in nature while the radiation associated 
with a beam of lower velocity is electrostatic in nature. 

For a very low beam velocity, (v/v ) cos 0 < 10, it is 

v XX 

difficult for the Cerenkov radiation to escape in a wide 
frequency range between the local plasma frequency and the 
local upper hybrid frequency due to Landau damping. Thus, 
the only escape process for waves of short wavelength 
(shorter than the limiting value for the Landau damping) 
corresponds to perpendicular mode electrostatic wave which 
propagates through the inhomogeneous plasma . 

CONCLUSIO N 

Numerical calculations were carried out for the disper- 
sion relation over a wide range of wavelengths from the 
electromagnetic waves to the electrostatic waves for a 
collisionless Maxwellian magnetoactive warm plasma. The 
computations were made mainly for the case of propagation 
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perpendicular to the superimposed magnetic field for 
several plasma conditions from rr/fj = 0.2 to 4.5. The 
results indicate that the dispersion curves in the short 
wave length range (electrostatic wave range) merge into 
the harmonics of the electron cyclotron frequency when 
the wave length becomes longer except for the case of the 
upper hybrid resonance frequency. The merging of the disper- 
sion curves into the harmonics of the electron-cyclotron 
frequency indicates that the electrostatic waves correspond- 
ing to these curves are subjected to strong cyclotron 
damping and thus cannot escape as electromagnetic waves. 

The upper hybrid frequency is neither the resonance frequency 
of the electromagnetic wave nor the cut-off frequency of 
the electrostatic wave, but designates the intermediate 
wave length region where the wave characteristics change 
gradually from electrostatic to electromagnetic with increas- 
ing wave length. When an electrostatic wave packet consisting 
of a narrow band of frequencies around the center frequency 
starts to propagate through a slightly inhomogeneous plasma 
in the direction of decreasing tt/Q, the wave packet is re- 
flected at the point where the group velocity changes from 
negative to positive, and propagates in the opposite direction. 
While passing through the upper hybrid resonance point, 
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the wave character gradually changes from 


I 2 2 

u) o 63 v 17 + ^ » 

electrostatic to electromagnetic. The special variation 

of the electron density near this point gradually changes 

the direction of the ray path. The wave normal direction 

also turns in the oblique direction with respect to the 

superposed magnetic field. Certain amounts of energy can 

be shared with the parallel component of the wave, i.e., 

the component with a wave normal directed parallel to the 

constant magnetic field. At the plasma frequency point 

(uj = tt) part of the energy can be divided into the O-mode 

and part of it into the X-mode. The O-mode escapes from 

the plasma, whereas the X-mode continues to propagate in 

the plasma until it is reflected at the oj s id cutoff 

o z 

point. When this reflected wave again encounters the = tt 
point, part of the wave energy escapes as an O-mode EM wave 
and part is converted to an electrostatic wave and is ab- 
sorbed by the plasma. 

2 2 

tt + Q > uj > Q (for Q > tt) or 

— 2 2 ~ 

y tt +0 > uj > tt (for tt > 0) , where the above described 

conversion of electrostatic wave energy to electromagnetic 
wave energy takes place, coincides with the frequency range 
of observed electromagnetic Cerenkov radiation. Thus this 
electromagnetic radiation may originate as electrostatic 
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Cerenkov radiation associated with a beam of low velocity 

electrons, (v,, < v < 100 v iU ) in addition to the electro- 

magnetic Cerenkov radiation associated with a high velocity 

beam (v > 100 v ±1 ) . It is also possible that plasma in- 

th 

stabilities give rise to electrostatic waves which are 
then converted to electromagnetic waves; these waves may 
be erroneously attributed to Cerenkov radiation. 
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FIGURE CAPTIONS 


Figure 1 0 

Figure 2. 
Figure 3. 
Figure 4. 
Figure 5. 


w/0 - kR diagrams (dispersion curves) of the 

plasma waves which propagate in the direction 

perpendicular to the constant magnetic field for 

q„ = 0.2 and 0.8; k is equal to k in the text. 

N x 

Same as Fig. 1 with q^ = 2.0 and 2.5. 

Same as Fig. 1 with q^ = 3.0 and 3.5. 

Same as Fig. 1 with q XT = 4.0 and 4.5. 

N 

Change of the relative position of the frequency 

range, corresponding to the test wave packet, 

from A , to E J while the wave packet progresses 
d a 

through an inhomogeneous plasma. 


Figure 6. Schematic illustration of the path of the wave 
packet propagation through the inhomogeneous 
plasma in which the gradient of plasma density 
vn makes an arbitrary angle with respect to the 
magnetic field B ; the points from A to E 
are discussed in the text corresponding to 
to E^ of Figure 5. 


j* 
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Figure 7. 


u)/Q - kR diagrams for propagation in the direction 


with angle 0° to 90° with respect to the magnetic 

-2 

field; Dotted area between kR = 9 x 10 to 

4 x 10 , where k = k /sin 0, indicates that the 

x 

waves propagating in oblique directions are sub- 
jected to damping. 

Figure 8. Graphical solution of equation 13 relating the 

wavelength of Cerenkov radiation to the particle 
velocity v; the solutions are indicated by the 
cross points between the solid and dashed curves 
corresponding to the given value of v and 0; 

0 dependency of the solid curves is same as 
Fig . 7 . 
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